Phantom evolving wormholes with big rip singularities 
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We investigate a family of inhomogeneous and anisotropic gravitational fields exhibiting a future 
singularity at a finite value of the proper time. The studied spherically symmetric spacetimes are 
asymptotically Friedmann-Robertson- Walker at spatial infinity and describe wormhole configura- 
tions filled with two matter components: one inhomogeneous and anisotropic fluid and another 
isotropic and homogeneously distributed fluid, characterized by the supernegative equation of state 
cj=p/p<— 1. In previously constructed wormholes, the notion of the phantom energy was used 
in a more extended sense than in cosmology, where the phantom energy is considered a homoge- 
neously distributed fluid. Specifically, for some static wormhole geometries the phantom matter was 
considered as an inhomogeneous and anisotropic fluid, with radial and lateral pressures satisfying 
the relations Pr/ p < — 1 and Pr, respectively. In this paper we construct phantom evolving 

wormhole models filled with an isotropic and homogeneous component, described by a barotropic or 
viscous phantom energy, and ending in a big rip singularity. In two of considered cases the equation 
of state parameter is constrained to be less than —1, while in the third model the finite-time future 
singularity may occur for a; < — 1, as well as for — 1 < a; < 1. 

PACS numbers: 04.20.Jb, 04.70.Dy,11.10.Kk 



I. INTRODUCTION 



Recent astrophysical observations indicate that our 
Universe is currently in accelerating expansion [l|, 
This discovery has stimulated an intensive study of mod- 
els where a large number of possible cosmological mech- 
anisms have been proposed to explain the origin of the 
current acceleration. In the framework of general rela- 
tivity, a quantitative analysis shows that a mysterious 
component of energy, dubbed dark energy, is respon- 
sible for the origin of this cosmic phenomenon and it 
starts dominating the matter content dynamics only at 
recent times, being irrelevant at earlier stages of the evo- 
lution. Many dark energy models have been proposed to 
solve this fundamental problem of cosmological physics, 
such as A-CDM model Q, quintessence models 0, 13] > 
k-essence models Q, phantom models @, etc. 

These new advances in cosmology allow us to consider 
new type of singularities, besides already considered stan- 
dard singularities such as Big-Bang and Big-Crunch Q. 
This non-standard type of singularities occurs at a finite 
value of cosmological time and are included in the fol- 
lowing descriptive classification: the singularity is a Big 
Rip when the scale factor, energy density and pressure 
go to infinity in a finite proper time, and is a sudden sin- 
gularity when at a finite value of time and scale factor. 
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curvature or one of its higher derivatives blow up [s"] . 

The main motivation of this type of singularities comes 
from phantom cosmological models. In the framework 
of Friedmann-RobertsonWalker (FRW) cosmologies filled 
with two matter contents and dominated by a phantom 
type fluid, satisfying a non-dissipative barotropic equa- 
tion of state, this new type of singularities can be com- 
pletely classified in the following four types: type I for 
a Big Rip singularity; type II for a sudden singularity 
defined by a finite energy density and diverging pressure; 
type III for diverging energy density and pressure at a 
finite value of the scale factor, and the type IV for finite 
curvature components and diverging higher derivatives 
oiH i. 

On the other hand, since the pioneering works by Mor- 
ris and Thorne T^ , the study of wormholes has become 
one of the most popular and intensively studied topics 
in relativistic physics, where most of the efforts are di- 
rected to study Lorentzian wormholes, in the framework 
of classical general relativity, sustained by an exotic mat- 
ter with negative energy density. These models include 
both static [ll] and evolving relativistic versions [l^, sus- 
tained by a single fluid component. The interest has been 
mainly devoted to traversable wormholes, which have no 
horizons, allowing two-way passage through them [l3j . 

For static wormholes the fluid requires the violation of 
the null energy condition (NEC), while in Einstein grav- 
ity there are nonstatic Lorentzian wormholes which do 
not require WEC violating matter to sustain them. Such 
wormholes niay exist for arbitrarily small or large inter- 
vals of time [14{ or even satisfy the dominant energy con- 
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dition (DEC) in the whole spacetime [H, [T^ . One can 
consider also dynamic wormhole spacetimes filled with 
two fluids, just like it is required in cosmology where 
such two-fluid models are widely considered today in or- 
der to explain the observed accelerated expansion of the 
Universe [l7| . 

Wormhole spacetimes filled with phantom type mat- 
ter were considered before [l^. Specifically, spheri- 
cally symmetric static wormholes were studied, sustained 
by a phantom type matter (or super-quintessence) with 
anisotropic pressure. In these static models the notion 
of the phantom energy is used in a more extended sense 
than in cosmology since, strictly speaking, the phantom 
matter is a homogeneously distributed fluid, and for these 
non-dynamic wormhole models an inhomogeneous and 
anisotropic matter component is used, having for the ra- 
dial and lateral pressures Pr < —1 and 7^ Pr, respec- 
tively [1^. This type of extended phantom-like matter 
was also used in dynamical wormhole models [20]. 

The time evolution of wormhole geometries in a Fried- 
mann universe exhibiting a Big Rip singularity were pre- 
viously studied. However the literature on this topic is 
not extensive. For example, in Ref. [2l| the author con- 
siders some accelerated higher-dimensional cosmologies 
with a traversable static wormhole, dominated by a time- 
dependent cosmological constant, and ending at a Big 
Rip. The studied Big Rip solutions have an exponential 
scale factor. In addition, the authors of Ref. [13 con- 
sider two different wormhole models, modeled by a thin 
spherical shell accreting the phantom fluid. 

It is interesting to note that in this context, it has re- 
cently proposed that it is possible that the universe could 
avoid the big rip singularity with the occurrence of a big 
trip [2^ [IJl , which is a cosmological event that may ap- 
pear during the evolution of a wormhole embedded in a 
FRW universe approaching the Big Rip singularity. In 
this case, the wormhole accreting phantom matter ex- 
pands faster than the background FRW universe, and 
the radius of the wormhole throat diverges before the 
Big Rip is reached. In this scenario, the wormhole en- 
gulfs the entire universe, which will reappear from the 
other wormhole throat |22| . 

In the present paper we intend to study evolving worm- 
holes filled with two matter components, where one of 
them is an isotropic homogeneously distributed phantom 
fiuid characterized the supernegative equation of state 
uj = p/p < —1, and presenting a future singularity at a 
finite value of the proper time. 

The organization of the paper is as follows: In Sec. 
II we present the dynamical field equations for worm- 
hole models with a matter source composed of an ideal 
isotropic cosmic fluid and an anisotropic and inhomoge- 
neous one. In Sec. Ill expanding wormholes filled with a 
barotropic dark energy and phantom fiuids are studied. 
We discuss explicit models ending at a finite-time future 
singularity. In Sec. IV viscous expanding wormholes are 
discussed. Models that may evolve to a finite-time future 
singularity are considered, and in Sec. V we conclude 



with some remarks. 



II. FIELD EQUATIONS 

In this paper we shall make use of some previously ob- 
tained results by one of the authors in Ref. [l^ . Let us 
state the main result obtained in that paper, concerning 
with solutions containing two fluids and admitting spher- 
ical symmetry in the framework of the Einstein gravity 
theory. 

Taking the metric 

ds' ^ -dt' + a{tf ^ ^ + r^dn^^ , (1) 

in comoving coordinates, filled with the anisotropic and 
inhomogeneous fluid p.^{t,r), and the isotropic and ho- 
mogeneous fluid p{t), the Einstein equations are given 

by 



orr? 3fc b' 

+ ^ 




= np^^{t,r) + np{t) + K,{2) 
= npr{t, r) + np{t) — A, (3) 
= Kp^{t,r) + Kp{t) - Kl'i) 



where dVi^ — dO"^ + sin?OdLp'^, k — SttG; A is the cos- 
mological constant, a{t) is the scale factor, k = —1,0, 1; 
H = a/a; and an overdot and a prime denote differenti- 
ation d/dt and d/dr respectively. 

In this case the 4- velocity of the fluids is given by the 
timelike vector u" = (1,0,0,0), and the radial and tan- 
gential pressures obey the barotropic state equations 

Pr{t,r) = UJrPi„{t,r), 

Piit^r) =u},p^^{t,r), (5) 

with constant state parameters ujr and oj, . 

Note that the essential characteristics of a wormhole 
geometry are encoded in the spacelike section of the met- 
ric (IlJ. It is clear that this metric becomes a zero-tidal 
force static wormhole if a{t) — >■ const, and as b{r) — > kr^ 
it becomes a fiat FRW metric for fc = 0, a closed FRW 
metric for k — 1, and an open FRW metric for k = —1. 

It can be shown that Eqs. ([I])-(|4]) may be rewritten in 
the form [3] 

ds^ = dt^ 

3i/2 + ^- = Kp(t)+A, (7) 
p + 3Hip + p) = 0, (8) 
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where C is an integration constant, and the inhomoge- 
neous and anisotropic cosmic fluid is given by 



Pr 

Pi 



(9) 
(10) 

(11) 



Here, the constraint ujr + 2a;, +1 = for Eqs. ([5]) was 
used. From these expressions we conclude that if w,, > 
or < — 1 the obtained gravitational configurations are 
asymptotically FRW solutions at spatial infinity. 

In conclusion, the main result of the Ref. is that 
the evolution of the scale factor a{t) in the metric © 
is governed by the standard Friedmann equations ([7]) 
and (O, and it is determined by the fluid p{t). This mat- 
ter component may be in principle an ideal barotropic 
fluid or any other cosmic fluid satisfying the requirements 
of isotropy and homogeneity. 



III. EXPANDING WORMHOLE UNIVERSES 
FILLED WITH A BAROTROPIC DARK 
ENERGY AND PHANTOM FLUIDS 

In this section we shall consider that the isotropic 
and homogeneous matter component is described by a 
barotropic phantom energy p{t) with an equation of state 
of the form 



pit) = 



(12) 



where the constant state parameter u satisfles the con- 
straint w < — 1. In general Eq. allows us to consider 
a barotropic matter component describing standard mat- 
ter for a; > 0, a dark energy fluid for — 1 < w < —1/3, 
and super-quintessence for w < — 1, just like it is defined 
in cosmology. 

From now on in this section, we shall consider solutions 
with fc = A = and k = SttG = 1. Thus in this case, 
from Eqs. ([7]) and ([8]), we obtain that the scale factor is 
given by a{t) = ^(F + (3/2)(w + l)^)^/^^^"+^^^ where D 
and F are constants of integration. This scale factor we 
shall rewrite as 



ait) 



1 



l)t 



2/(3(c.. + l)) 



and the energy density in the form 

Po 



Pit) 



1 



^o(^ + lK 



2 ' 



(13) 



(14) 



where — 3iJ^, in order to have a{tg = 0) = a„ > 
and Hit^ = 0) = iJ„ > 0. 

It is easy to verify that in this case the metric (|6]) 
becomes 

4/(3(<^ + l)) 



ds^ = - aM 1 + -H^ (w + l)t 



1 - 



i^) 



(15) 



and the anisotropic and inhomogeneous energy density is 
given by 



P.At,r) 



(^) 



-(l+3a;r)/LJr 



H„{^ + l)t 



(16) 



The metric (jl5|) represents an evolving wormhole with a 
throat located at tq for < — 1 and ojr > fl6l] and is 
asymptotically a flat FRW universe. 

It is well known that in general to keep a wormhole 
open exotic matter with a negative energy density at 
the throat is needed [l3|. However, there are examples 
of evolvingwormholes satisfying the DEC in the whole 
spacetime jl5l [l6j , which implies that the energy density 
is positive everywhere. In the studied here solutions the 
branch with > is characterized by a positive radial 
pressure pr, and negative p.^ and pi. But, there are also 
wormholes which avoid the usual exotic matter require- 
ments for wormholes. The branch with Wr- < — 1 has a 
positive energy density and lateral pressure, while the 
radial pressure is negative and larger in magnitude than 
the energy density. Speciflcally, the total matter content 



p^(i,r) = p(<) + p,„(t,r) 



(17) 



for these models is determined by Eqs. and p^ . 
Thus, for any value of the state parameters lo and < 
we have that always /o^ > 0. For ujr > we can have in 
general time intervals where the total energy is positive 
or negative. In effect, for uj > —1/3 the wormhole model 
starts with a positive total energy density (since for a 
flxed value r = const the isotropic component dominates 
over the another one), then decreases till zero at certain 
to, and becomes negative for t > to. For cj < —1/3 the 
total energy density starts negative, then increases till 
zero at certain to, and becomes positive for t > to- If 
Lj — —1/3 both components of total energy density ([T7]) 
behave as ~ 1/a^, and the expansion occurs at a constant 
velocity. 

In order to have accelerated (decelerated) expansion 
we must consider a;<— l/3(w> —1/3). Notice that at 
t = we have that the total energy density is given by 



-(l+3l.Jr)/tJr 



(18) 



We can see from Eq. (|T8| that if the values of Ur are 
constrained to be in the ranges w,. < —1/3 and iOr > 
0, then the relation p{0) > p.^ (0, r) is fulfilled for r > 



i- 



'po'Y 



not fulfilled, then at t 



r/iWujr)^ If the last inequality is 
the inhomogeneous component 
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FIG. 1: Penrose diagram for an evolving wormhole with uur < 
— 1 and uj < —1. The dotted line is the wormhole throat and 
the dashed line represents the future Big Rip singularity. 



begins dominating. If —1/3 < ujr < we have that 



ft„(0,r) > p(0) for r > {-rl+^^^'-uj-^a-^ p^^r-/(^+^'^-\ 
For these models it is useful to state that if the isotropic 
perfect fluid satisfies the DEC, i.e. — 1 < w < 1, we 
can rescale the cosniological time so that 1 + |_ffj,(a; + 
l)t — > t and the solution for the scale factor takes the 
standard form a{t) = aot'^^^-^'--^~^^^\ with the isotropic 
energy density given by p{t) — po/t^. 

In order to study these models during a phantom evo- 
lution we must consider the values w < — 1. In this case 
we have that l + ^Hg{uj + l)t = 1— ^H^\uj + l\t, so during 
the cosmic evolution this expression may vanish at some 
value of time t > 0. Thus we conclude that if w < — 1 we 
have a future singularity at a finite value of the proper 
time 



(19) 



since a{t) — > cxd, p{t) — >■ oo and p — > — oo at < = t^^^. 
By considering only the behavior of scale factor and the 
isotropic fiuid we may conclude that this future singular- 
ity is of a Big Rip type. 

However we must also consider the presence of the in- 
homogeneous matter component (|16p . As we have stated 
above in this case we have that p.^ (t, r) > for Ur < 0, 
and p.^{t^r) < for uJr > 0. For uJr — ~l/3 this matter 
component becomes homogeneous and isotropic one with 
~ ~ ~p.^j3. Note that for w > —1, uJr < —1/3 
and > we have that p^^(tc,r — > oo) — > and 
Pini^ — ^ oOjTc) — > for constant values tc and Vc- 

For Lj < —1, UJr < —1/3 and > we have that at 
t = the anisotropic and inhomogeneous energy density 
starts to evolve from the constant value p.^{0,rc), where 
the constant Tc > rp, and then decreases, becoming zero 
at t = t^^ . 

As a consequence, for t = t^^ and u < —1, we have 
that a{t) — !• oo, p{t) — !• oo, p — > — oo, p.^^(t,r) = 0, 
p^(t, r) — (t, r) — 0. This allows us to conclude that if 



we consider the total matter content given by Eq. pTl) . 
for t = t^^ and u < —1, we have that a{t) — > oo, 
PT{t,r) — !■ oo, PT^{t,r) — > -oo, PTi{t,r) — > -oo, 
where p^^{t,r) — p{t) + p^{t,r) and Pr^i{t,r) = pit) + 
Pi{t,r). In conclusion, the scale factor, the total energy 
density and the total pressures blow up at the finite time 
tj,^, so this finite-time future singularity is of a Big Rip 
type. In Fig.[T]is shown a conformal diagram of an evolv- 
ing wormhole with < — 1 and a phantom energy with 
w < -1. 

Lastly, let us introduce the equilibrium time defined by 
the condition p(t^^ ) — p.^^ (t^^ , r) , for wormholes models 
with a positive total energy density and big rip. This 
allows us to find i „ as a function of the radial coordinate 



r: 



3F„(l+a;) 



ir \ -2TT 



V 



r'^uJralpo 



J 



.(20) 



It becomes clear that in general teq can take complex 
as well as real values. In order to have always real values 
of the equilibrium time the condition a;^ < must be 
required. For wormhole models with a positive total en- 
ergy density and big rip the equilibrium time is positive 
ifa;r<— l,w<— 1 and 



(S) 



> 1. 



Thus, for To < r < r*, where 



(21) 



(22) 



we have that t^^ > 0, and t^^ < for r > r*. No- 
tice that from Eqs. ((T8| and ((2T|) we obtain for the total 
energy density the equivalent constraint Py(0,r) > 2po, 
and from Eqs. and (PO)) the condition t^q < ti,r is 
automatically fulfilled for < 0. By using Eqs. (PO]) - 
((22)) qualitative plots for wormholes defined by conditions 
< — 1 and (jjr > are shown in Figs. 2(a)|2(c)" and in 



Figs.[3(iOp(^ 



respectively. 

Now we want to study evolving wormhole models with 
future singularities and filled with a viscous phantom 
matter. This issue will be addressed in the next section. 



IV. VISCOUS EVOLVING WORMHOLES 

Let us now consider wormhole models with the 
isotropic and homogeneous matter component described 
by a viscous phantom fluid. The role of the dissipa- 
tive processes has been extensively considered in cosmol- 
ogy |25l . l26j , where the study is done within the frame- 
work of the standard Eckart theory of relativistic irre- 
versible thermodynamics. Any dissipation process in a 
FRW cosmology is scalar, and therefore may be modeled 
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(a) ro < Tc < r» 



(b) rc > r. 



(c) Tc > ro (r* < ro) 



FIG. 2: These figures show, for wormhole models with Ur < —1 and iu < —1, the qualitative behavior of the dimensionless 
energy density for the isotropic component Q,{t) (dotted line), anisotropic matter i^i„{t,rc) (dash-dotted line) and the total 
matter content nT{t,rc) (solid line) as a function of the time t, at a constant value of the radial coordinate Vc- The comoving 
time varies from to the big rip time f^^, . In Figs. 2(a) and 2(b) are plotted curves fulfilling the constraint (I21|l : if ro < < r. 



the anisotropic component dominates for < t < t^^, while the isotropic component dominates over the anisotropic one at 
the time interval t^^ < t < t^^ (Fig. 2(a) I; on the other hand if rc > r* the isotropic component always dominates over the 
anisotropic for < t < (Fig- |2(b)[ ). In Fig. |2(c)| the constraint (|2ip is not fulfilled, thus always the isotropic component 
dominates over the anisotropic one for Q < t < t^^ and any constant value r-c > ro. 



n(t) mii(t,r) nt(t,r) Q(t) mnltj) Qt(t,r) Q(t) ilinltj) Qt(t,r) 



a 




FIG. 3: These figures show, for wormhole models with uj,- > and uj < —1, the qualitative behavior of the dimensionless energy 
density for the isotropic component r2(t) (dotted line), anisotropic matter ^in{t,rc) (dash-dotted line) and the total matter 
content Q,T{t,rc) (solid line) as a function of the time t, at a constant value of the radial coordinate rc. The comoving time 
varies from to the big rip time t^^. In this case always the isotropic component dominates over the anisotropic one. The 
total energy density vanishes at t = Iq. For Fig. 3(a) we have taken < to < tbr, while for Figs 3(b) and 3(c) to < 0. Here 



: ro(roajraopo) 



Note that in this case the anisotropic energy density is always negative. 



as a bulk viscosity within a thermodynamical approach. 
The bulk viscosity introduces dissipation by only redefin- 
ing the effective pressure, P.^.^. , according to 

P^j^ = p -I- n = p - ^H, (23) 

where 11 = Ii{t) is the bulk viscous pressure, ^ — ^{t) is 
the bulk viscosity coefficient and H is the Hubble param- 
eter. 

In this case the Friedmann equations (O and ([S]), with 



k — and k — SirG — 1, take the form 

ZH^ = p + A, (24) 
p + 3H(p + p + n) = 0. (25) 

The violation of DEC is expressed by the relation p~\-p + 
n < 0. This condition implies an increasing energy den- 
sity of the isotropic fluid filling the evolving wormhole, 
for a positive bulk viscosity coefficient. The condition 
^ > guaranties a positive entropy production and, in 
consequence, no violation of the second law of the ther- 
modynamics [27| . 
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We shall assume that the viscous component obeys the 
state equation ([T^ . hence from Eq. (1^^ we have that 
P^jj = up - 3^H. Thus from Eqs. 1^ and ^ we 
obtain the following evolution equation for H: 

2H + 3{uj + l)H^ = 3^H +{u + l)A. (26) 

From this equation we obtain for A = that 



Hit) 



Thus for 7^ — 1 the scale factor is given by 

/ o r , r \ 2/(3(1^ + 1)) 

while for a; = — 1 it may be written as 



(27) 



i{t) = De- 



dt 



(28) 



(29) 



where C and D are integration constants. In general, for 
the considered case, the solution may be written through 
^(t) or a{t) because there are three independent equations 
for the four unknown functions a{t), p(t), ^(t) and p{t). 
In our case we have written the solution through the bulk 
viscosity £,{t). It is worth to mention that for a given a{t) 
we can write H and then obtain the expressions for the 
energy density from Eq. (j24p and the bulk viscosity from 

Eq. dini). 

It becomes clear that for ^ = and w 7^ — 1 we obtain 
from Eq. (l28l) the solution (fT3)) discussed in the previous 
section, while for a vanishing bulk viscosity the de Sitter 
scale factor a{t) — e^o* is obtained for p — —p — const. 

On the other hand, note that from Eq. and 
Eq. (j25p we may write that 



^^h + h^^-Up- 

a 6 



3P 



A 
3"' 



(30) 



Thus the condition for an expansion with constant ve- 
locity is given by p-|-3P^^.^ = 2A. By taking into account 
Eq. (j23p we may write 



^((l + 3.)p-2A) 



(31) 



Note that for A = 0, and by taking into account Eq. ([M]). 
we conclude that in order to have dynamic wormholes 
expanding with constant velocity the bulk viscosity must 
be given by 



(1 



3\/3 



3w) 1/2 
P ■ 



(32) 



In this case we see that a necessary condition to have a 
positive bulk viscosity coefficient is that uj > —1/3. 

Now we shall consider specific viscous phantom evolv- 
ing wormhole models. 



A. Wormhole models with constant bulk viscosity 

Let us now consider wormhole models with a vanishing 
cosmological constant and a bulk viscosity given by 



^{t) = Co — const. 



(33) 



For a; 7^ —1, the Eq. ((28l) allows us to write the scale 
factor in the form 



2/(3(1^ + 1)) 



(34) 



from which we obtain that the homogeneous and 
isotropic energy density is given by 



Pit) = 



(^l + |i(tj + l)(e3?ot/2_l)y 



(35) 



Then the anisotropic and inhomogeneous matter compo- 
nent takes the following form: 



'«Pi„ (*> 



-(l+3Wr)/tJr 



2c.,a2(l + ^(w-)-l)(e3?o */2_i)) 



4/(3(<.. + l)) 



(36) 



As in the previous section, for a; < — 1 we have a future 
singularity at a finite value of the comoving proper time 
tj_^ since a(t) — > 00, p{t) — > 00 and p — > —00 at 



t 



«0 



> 0. In this case at t—t. the 



energy density of the anisotropic matter threading the 
wormhole vanishes since if > or < ^ 1 we obtain 
that p.^{t^^,r) — for any r > tq. 

If we consider the total energy density , given now 
by Eqs. ^ and we conclude that for w < -1 
— > 00, Pj,^ — > 00 and — > 00 at t = t^^. Thus, 
this future singularity is characterized by diverging scale 
factor, total energy density and total pressures, but with 
a well behaved bulk viscosity, since ^ is constant during 
all evolution. 

Notice that for the lj — —1 branch solution we obtain 
from Eq. that the scale factor is given by 



a{t) = De 3«o 



(37) 



In this case the model is characterized by an accelerated 
expansion and does not end in a future singularity. 



B. Accelerating wormhole models with ^ ~ p^^^ 

Another interesting example in this line is obtained 
for a bulk viscosity given by ^ = ap^^'^, where a is a 
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constant parameter. In this case, for any value of the 
state parameter lu and A = 0, the integration of Eq. ([25)1 
allows us to write 



H 



Ho 



l + |i?o(t^ + l-\/3a)i' 



(38) 



where Hq = H{t = 0). Thus, the scale factor becomes 



a big rip singularity not only for viscous phantom energy, 
but also for viscous dark energy, and even for standard 
viscous matter (see Figs. 4(a) and |4(b)[ ). Clearly all these 
models have an accelerated expansion. It is easy to show 
that in order to have models expanding with constant 
velocity we must require that — — = 1. This im- 

plies that a = {1 + 3a;)/3\/3 in agreement with Eq. (|32p . 



a{t) = ao 1 + 7:Hq{u + 1 - V3a)t , (39) 



with flo = a{t = 0). The energy density of the isotropic 
component takes the form 



p{t) = 



3iJ2 



(l + |i7o(t^ + l-\/3a)t)2' 



(40) 



while the bulk viscosity and the energy density of the 
anisotropic and inhomogeneous fluid are given by 



1 + ^Ho{io + 1 - V3a)t' 



(41) 



-(l+3wr)/i.Jr 



r2w^a2 (1 + ^Ho{uj + 1 - Via)t) 3(-+i-v^.) 



(42) 



respectively. 

By demanding that 



VSa > a; + 1 



(43) 



the scale factor, the isotropic energy density and pressure 
blow up to infinity at a finite time 



2H' 



br 



3(V3a - (tJ + 1)) 



>0, 



(44) 



and then we have the occurrence of a future singularity. 
As in the previous case, at t=ij^^ the energy density and 
pressures of the anisotropic matter threading the worm- 
hole vanishes for any r > tq. 

By considering the total energy density , given now 
by Eqs. ^ and (gH), we conclude that for the 

constraint (j43|) we have that — > oo, p^^ — > oo and 
Py, — > oo at t = t^^. Thus, this future singularity is 
characterized by diverging scale factor, total energy den- 
sity and total pressures, and diverging too bulk viscos- 
ity (1411) . since it blows up at the time if,,-. In conclusion, 
in these viscous expanding wormhole models the scale 
factor, the total energy density, the total pressures and 
the bulk viscosity blow up at the finite time t^^, so this 
finite-time future singularity is of a Big Rip type. 

It is interesting to note that the constraint P5)) implies 
that if the bulk viscosity is positive, i.e. a > 0, we can 
have a future singularity also for uj > —1. Thus we have 



V. CONCLUSIONS AND FURTHER 
COMMENTS 

This paper deals with inhomogeneous and anisotropic 
spacetimes, ending in a future singularity at a finite value 
of the proper time, and filled with an inhomogeneous 
and anisotropic fiuid and another isotropic and homoge- 
neously distributed super-quintessence fluid. The studied 
solutions describe evolving wormholes for which the rate 
of expansion is determined by the phantom energy, while 
the inhomogeneous and anisotropic component threads 
and sustains the wormhole. The main purpose of this 
work is to present analytic wormhole models exhibiting 
a Big Rip during its evolution. Three independent cos- 
mological models are explored. 

In the first model it turns out that, the isotropic and 
homogeneous component is a barotropic phantom fiuid. 
For these evolving wormholes the scale factor, the total 
energy density and the total pressures blow up at a finite 
proper time, so this finite-time future singularity is of a 
Big Rip type. In the second wormhole configuration the 
isotropic and homogeneous component is a viscous phan- 
tom fiuid with a constant bulk viscosity. Now, the future 
singularity is characterized by diverging scale factor, to- 
tal energy density and total pressures, but with a well 
behaved bulk viscosity, due to the constant character of 
^ during all evolution. In the latter model the isotropic 
and homogeneous component is a viscous phantom fluid 
with a bulk viscosity of the form ^ ~ p^^"^. For these 
viscous dynamic wormholes the scale factor, the total 
energy density, the total pressures and the bulk viscosity 
blow up at a finite proper time, so this future singularity 
is also of a Big Rip type. It must be added that this 
third model allows us to consider Big Rip wormholes not 
only for a viscous super-quintessence energy. Effectively, 
if ^ > 0, we can have a Big Rip also for viscous dark 
energy (i.e. for —l<uj< —1/3), and even for standard 
viscous matter (i.e. for a; > —1/3). If a < the Big 
Rip is avoided. However in this case the bulk viscosity 
becomes negative and, in consequence, the second law of 
the thermodynamics is not fulfilled. 

Notice that, in all considered here solutions the mixed 
component of the energy-momentum tensor Ttr vanishes. 
This means that there is no radial energy flow and no 
accretion onto the wormhole of phantom energy from 
the cosmic fluid. Thus the mechanism by which the big 
trip could be achieved is out of the possibilities for these 
wormhole models 1241. 



8 




for the isotropic component fl{t) (dotted line), anisotropic matter fli„{t,rc) (dash-dotted hne) and the total matter content 
QT{t,rc) (solid line) as a function of the time t, at a constant value of the radial coordinate rc. The comoving time varies 
from to the big rip time t^^. For wormholes defined hy Ur < ~1, we can have a similar behavior for all values of the state 
parameter uj, i.e. for oj < —1, — 1 < a; < —1/3 and uj > —1/3, as is shown in Fig. 4(a) 
^^in(i, rc) ~ Q{t) and the big rip singularity is present even if cj > —1 

behavior for all values of the state parameter u, i.e. for uj < —1, —1 < u < —1/3 and uj > 

that in this case the anisotropic energy density is always negative. The big rip singularity is also present even if uj > 
to represents the time at which QT{t,rc) = 0. 



Here t„ 



represents the time at which 
For wormholes defined by cJr > we have again a similar 

Note 



-1/3, as is shown in Fig. 4(b) 



-1 and 



Lastly, we want to state that all obtained here results 
on future singularities are applicable also for flat FRW 
cosmological models, since all the discussed wormhole so- 
lutions are asymptotically flat FRW cosmologies. 



tigacion de la Universidad del Bio-Bi'o through grants N° 
DIUBB 121007 2/R and GI121407/VBC (MC). 



VI. ACKNOWLEDGEMENTS 



This work was supported by CONICYT through Grant 
FONDECYT 1080530 and by the Direccion de Inves- 



[1] D. A. Howell, M. Sullivan, P. E. Nugent, R. S. EUis, 
A. J. Conley, D. Le Borgne, R. G. Carlberg and J. Guy 
et al, Nature 443, 308 (2006); A. G. Riess et al. [Su- 
pernova Search Team Collaboration], Astrophys. J. 607, 
665 (2004); J. L. Tonry et al. [Supernova Search Team 
CoUaboration], Astrophys. J. 594, 1 (2003); S. Perlmut- 
ter et al. [Supernova Cosmology Project Collaboration], 
Astrophys. J. 517, 565 (1999); A. G. Riess et al. [Super- 
nova Search Team Collaboration], Astron. J. 116, 1009 
(1998). 

[2] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. 
Phys. D 15, 1753 (2006). 

[3] A. V. Astashenok, S. 'i. Nojiri, S. D. Odintsov and 
R. J. Scherrer, Phys. Lett. B 713, 145 (2012); R. A. Suss- 
man and G. Izquierdo, Class. Quant. Grav. 28, 045006 
(2011); R. -J. Yang and S. N. Zhang, Mon. Not. Roy. 
Astron. Soc. 407, 1835 (2010); A. Kurek and M. Szyd- 
lowski, Astrophys. J. 675, 1 (2008); V. Sahni, Lect. Notes 
Phys. 653, 141 (2004); M. S. Turner and M. J. White, 
Phys. Rev. D 56, 4439 (1997). 



[4] G. Olivares, F. Atrio-Barandela and D. Pavon, Phys. 
Rev. D 71, 063523 (2005); V. F. Cardone, A. Troisi and 
S. Capozziello, Phys. Rev. D 72, 043501 (2005); T. Matos 
and L. A. Urena-Lopez, Phys. Rev. D 63, 063506 (2001); 
E. J. Copeland, N. J. Nunes and F. Rosati, Phys. Rev. 
D 62, 123503 (2000); P. F. Gonzalez-Diaz, Phys. Rev. D 
62, 023513 (2000). 

[5] M. U. Farooq, M. Jamil and M. A. Rashid, Int. J. Theor. 
Phys. 49, 2278 (2010); N. Cruz, P. F. Gonzalez-Diaz, 
A. Rozas-Fernandez and G. Sanchez, Phys. Lett. B 679, 
293 (2009); L. P. Chimento, Phys. Rev. D 69, 123517 
(2004); T. Chiba, Phys. Rev. D 66, 063514 (2002). 

[6] K. Karami, S. Ghaffari and J. Fehri, Eur. Phys. J. C 64, 
85 (2009); H. Mohseni Sadjadi and M. Alimohammadi, 
Phys. Rev. D 74, 043506 (2006); J. Kujat, R. J. Scher- 
rer and A. A. Sen, Phys. Rev. D 74, 083501 (2006); 
L. Perivolaropoulos, Phys. Rev. D 71, 063503 (2005). 

[7] O. -C. Stoica, Commun. Theor. Phys. 58, 613 (2012); 
G. Niz and N. Turok, Phys. Rev. D 75, 026001 (2007); 
C. Cattoen and M. Visser, Class. Quant. Grav. 22, 4913 



9 



(2005); J. M. M. Senovilla, Phys. Rev. Lett. 64, 2219 
(1990); G. L. Murphy, Phys. Rev. D 8, 4231 (1973). 

[8] P. Singh, Class. Quant. Grav. 26, 125005 (2009); 
M. P. Dabrowski and A. Balcerzak, gr-qc/0701056 

[9] S. Nojiri, S.D. Odintsov and S. Tsujikawa, Phys. Rev. D 
71, 063004 (2005). 
[10] M.S. Morris and K.S. Thorne, Am. J. Phys. 56, 395 
(1988); M.S. Morris, K.S. Thorne and U. Yurtsever, 
Phys. Rev. Lett. 61, 1446 (1988); M. Visser, Lorentzian 
Wormholes: From Einstein to Hawking, (AIP, New York, 
1995). 

[11] A. Anabalon and A. Cisterna, Phys. Rev. D 85, 084035 
(2012); G. A. S. Bias and J. P. S. Lemos, Phys. Rev. D 

82, 084023 (2010); A. B. Balakin, J. P. S. Lemos and 
A. E. Zayats, Phys. Rev. D 81, 084015 (2010); M. Jamil, 
P. K. F. Kuhfittig, F. Rahaman and S. A. Rakib, Eur. 
Phys. J. C 67, 513 (2010); M. Jamil and M. U. Farooq, 
Int. J. Theor. Phys. 49, 835 (2010); E. F. Eiroa, Phys. 
Rev. D 80, 044033 (2009); M. Cataldo, P. Salgado and 
P. Minning, Phys. Rev. D 66, 124008 (2002). 

[12] M. Cataldo, P. Meza and P. Minning, Phys. Rev. D 

83, 044050 (2011); S. H. Hendi, J. Math. Phys. 52, 
042502 (2011); S. H. Hendi, Can. J. Phys. 89, 281 

(2011) ; I. Bochicchio and V. Faraoni, Phys. Rev. D 82, 
044040 (2010); M. Jamil and M. Akbar, arXiv:0911.2556 
[hep-th]; M. Jamil, arXiv:0909.3689 [gr-qc]; J. Hansen, 
D. i. Hwang and D. h. Yeom, JHEP 0911, 016 (2009); 
M. H. Dehghani and S. H. Hendi, Gen. Rel. Grav. 
41, 1853 (2009); A. DeBenedictis, R. Garattini and 
F. S. N. Lobo, Phys. Rev. D 78, 104003 (2008); 
P. F. Gonzalez-Diaz, Phys. Rev. D 68, 084016 (2003). 

[13] F. S. N. Lobo and M. Visser, Class. Quant. Grav. 21, 
5871 (2004). 

[14] S. Kar, Phys. Rev D 49, 862 (1994); S. Kar and D. 

Sahdev, Phys. Rev D 53, 722 (1996). 
[15] H. Maeda, T. Harada and B.J. Carr, Phys. Rev. D 79, 

044034 (2009). 

[16] M. Cataldo and S. del Campo, Phys. Rev. D 85, 104010 

(2012) . 

[17] M. Cataldo, P. Mella, P. Minning and J. Saavedra, Phys. 
Lett. B 662, 314 (2008); A. Banijamah, B. Fazlpour, 
Phys. Lett. B703, 366-369 (2011); S. Chattopadhyay, 
U. Debnath, Astrophys. Space Sci. 326, 155-158 (2010); 
H. FarajoUahi, N. Mohamadi, H. Amiri, Mod. Phys. Lett. 
A25, 2579-2589 (2010); Y. -F. Cai, J. Wang, Class. 



Quant. Grav. 25, 165014 (2008); L H. Brevik, Gen. Rel. 
Grav. 38, 1317-1328 (2006). 

[18] M. Jamil, P. K. F. Kuhfittig, F. Rahaman and 
S. .ARakib, Eur. Phys. J. C 67, 513 (2010); M. Jamil 
and M. U. Farooq, Int. J. Theor. Phys. 49, 835 (2010); 
J. A. Gonzalez, F. S. Guzman, N. Montelongo-Garcia and 
T. Zannias, Phys. Rev. D 79, 064027 (2009); M. Jamil, 
U. Farooq and M. A. Rashid, Eur. Phys. J. C 59, 
907 (2009); K. A. Bronnikov and A. A. Starobinsky, 
JETP Lett. 85, 1 (2007); F. Rahaman, M. Kalam, 
B. C. Bhui and S. Chakraborty, Phys. Scripta 76, 56 
(2007); P. K. F. Kuhfittig, Class. Quant. Grav. 23, 5853 
(2006); O. B. Zaslavskii, Phys. Rev. D 72, 061303 (2005). 

[19] F. S. N. Lobo, Phys. Rev. D 71, 084011 (2005); 
S. V. Sushkov, Phys. Rev. D 71, 043520 (2005). 

[20] M. Cataldo, P. Labrafia, S. del Campo, J. Crisostomo and 
P. Salgado, Phys. Rev. D 78, 104006 (2008); M. Cataldo, 
S. del Campo, P. Minning and P. Salgado, Phys. Rev. D 
79, 024005 (2009). 

[21] A. R. El-Nabulsi, Gen. Rel. Grav. 43, 261 (2011). 

[22] V. Faraoni and W. Israel, Phys. Rev. D 71, 064017 
(2005). 

[23] A. V. Astashenok, S. Nojiri, S. D. Odintsov and 
A. V. Yurov, Phys. Lett. B 709, 396 (2012); A. V. Yurov, 
A. V. Astashenok and V. A. Yurov, Astrophys. Space 
Sci. 342, 1 (2012); P. F. Gonzalez-Diaz, Phys. Lett. B 
635, 1 (2006); J. A. Jimenez Madrid, Phys. Lett. B 
634, 106 (2006); A. V. Yurov, P. Martin Moruno and 
P. F. Gonzalez-Diaz, Nucl. Phys. B 759, 320 (2006); 
P. F. Gonzalez-Diaz, Phys. Lett. B 632, 159 (2006); 
P. F. Gonzalez-Diaz and J. A. Jimenez-Madrid, Phys. 
Lett. B 596, 16 (2004); P. F. Gonzalez-Diaz, Phys. Rev. 
Lett. 93, 071301 (2004) 

[24] V. Faraoni, Phys. Lett. B 647, 309 (2007). 

[25] R. Colistete, J. C. Fabris, J. Tossa and W. Zim- 
dahl, Phys. Rev. D 76, 103516 (2007); A. Pradhan 
and P. Pandey, Astrophys. Space Sci. 301, 127 (2006); 
I. H. Brevik, O. Gorbunova and Y. A. Shaido, Int. J. 
Mod. Phys. D 14, 1899 (2005); I. H. Brevik and O. Gor- 
bunova, Gen. Rel. Grav. 37, 2039 (2005); W. Zimdahl, 
Phys. Rev. D 53, 5483 (1996). 

[26] M. Cataldo, N. Cruz and S. Lepe, Phys. Lett. B 619, 5 
(2005). 

[27] W. Zimdahl and D. Pavon, Phys. Rev. 61, 108301 (2000). 



